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Abstract: 

I consider a general block-tridiagonal matrix and the corresponding transfer 
matrix. By allowing for a complex Bloch parameter in the boundary conditions, 
the two matrices are related by a spectral duality. 

As a consequence, I derive some analytic properties of the exponents of the 
transfer matrix in terms of the eigenvalues of the (non-Hcrmitian) block matrix. 
Some of them are the single-matrix analogue of results holding for Lyapunov 
exponents of an ensemble of block matrices, that occur in models of transport. 
The counting function of exponents is related to winding numbers of eigenval- 
ues. I discuss some implications of duality on the distribution (real bands and 
complex arcs) and the dynamics of eigenvalues. 

P.A.C.S.: 02.10.Yn (matrix theory), 72.15.Rn (localization effects), 72.20.Ee 
(mobility edges, hopping transport) 



1. Introduction 

A tridiagonal Hermitian matrix whose entries are square matrices of size M 
is a block tridiagonal matrix. By denoting Hi = hJ the diagonal blocks and by 
Li the blocks in the next upper diagonal, the eigenvalue equation written in the 
block components of an eigenvector u is: 

H,Ui + Uui+i + = Eu, (1) 

We may view Hi as Hamiltonian matrices of a chain of subsystems, each with 
M internal states, sequentially coupled by matrices Li. We shall make the only 
requirement det Li ^ 0. 

Band matrices are in this class, with non-diagonal blocks being triangular. Ran- 
dom band matrices are studied for quantum chaos and transport, mainly by nu- 
merical means [1] or by a mapping on a non-linear supersymmetric sigma model 
[2]. The block structure is typical of tight binding models, as Anderson's model 
for the transport of a particle in a lattice with impurities [3] . Here, the matrices 
Hi are the Hamiltonians of isolated slices transverse to some direction and the 
matrices Li contain the hopping amplitudes between lattice sites of neighboring 
slices. Tridiagonal arrays of random matrices were studied in the context of 
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multichannel scattering in mesoscopic or nuclear physics [4] or as multi-matrix 
models in the large size limit [5]. For recent reviews of applications of random 
matrices in physics, sec [6, 7]. 

The second-order recursive character of (1) makes it useful to introduce a 
transfer matrix Ti{E), of size 2M x 2M 



J \u. 



(2) 



The block components of the eigenvector are reconstructed by applying a prod- 
uct of transfer matrices to an initial block pair. A length N of the chain corre- 
sponds to the transfer matrix T{E) =Tn{E)... Ti{E): 

The transfer matrix is the main tool for investigating boundary properties of 
the Hamiltonian's eigenvectors, or the transmission matrix of the chain coupled 
to leads. I developed a theory for transfer matrices of general block-tridiagonal 
Hamiltonians in [8, 9]. Recently the formalism was applied to transport in 
nanotubes and molecules and generalized to allow for non-invertible off-diagonal 
blocks [10, 11, 12, 13]. 

In this paper I describe some interesting consequences of a nice algebraic identity 
involving the characteristic polynomials of the two matrices [8]. Though this 
spectral duality holds in general, I here restrict to Hermitian block tridiagonal 
matrices because of their relevance in physics. 

For a chain of length we must provide boundary conditions. With periodic 
boundary conditions we require Lq = Ln for Hermiticity and Uq — U]^, Hn+i = 
Ui- However, it turns to be very convenient to allow for a complex Bloch 
parameter: 

Uq = ^Wjv. Hn+i = zui (4) 
Therefore the Hamiltonian matrix is block-tridiagonal with corners: 



n{z) 



2 

4 



Ln-1 

\zLn i]v-i / 



(5) 



It is Hermitian only for \z\ = 1. Boundary conditions with z = e^'^ arise when 
decomposing the eigen-problem (1) for an infinite periodic chain of period N in 
the eigen-spaces of the A-block translation operator, as well as in the topology 
of a A-site ring with a magnetic flux through it. In general it is: 

n{z)^='H{l/z*) (6) 
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Non Hermitian tridiagonal matrices (M = 1), with z = e^^ and <? > 0, were 
introduced by Hatano and Nelson [14] in a study on vortex depinning in super- 
conductors which promoted a burst of research, sec for example [15, 16, 17, 18]. 
Spectral properties were analyzed in greater detail by Goldsheid and Kho- 
ruzhenko [19] who proved that, for g > Qcr, eigenvalues corresponding to ex- 
tended states begin to migrate in the complex plane and distribute along the 
level curve of the single Lyapunov exponent of the model: j{E) = g. Two wings 
of real eigenvalues correspond to exponentially localized eigenvectors, that are 
insensitive to the boundary. 

These features also appear in the more difficult case of block matrices. For 
1^1 sufficiently greater or smaller than unity, the block matrix develops complex 
eigenvalues which arc seen numerically to distribute along lines. The eigenval- 
ues of an eptadiagonal matrix (M=3) with diagonal disorder and unit hopping 
amplitudes, are shown in Fig. 2. 

It is intuitive that there is a connection between boundary properties of 
eigenvectors, which are described by the transfer matrix, and the response of 
energy eigenvalues to variations of boundary conditions: namely, the Landauer 
and Thouless approaches to transport. This intuition has a formulation in the 
spectral duality. 

After a short review of duality (section 2), I derive several analytic properties 
of the exponents of a single transfer matrix (section 3). In particular, I evaluate 
the counting function of exponents as the winding number of trajectories of 
eigenvalues of the source block matrix. In section 4, I use duality to describe 
qualitatively the distribution and dynamics of eigenvalues of Hermitian and 
non-Hermitian block-tridiagonal matrices. 



2. Spectral Duality 

In this section I review some basic properties of a transfer matrix T{E) of a 
chain of length N, with Lq = Ln [8, 9]. 

T{E) is a matrix polynomial in E of degree TV, with a non-zero determinant 
independent of E and of diagonal blocks 11^: 



The block structure of the Hamiltonian and the corresponding factorization of 
the transfer matrix, makes a typical property of transfer matrices apparent: 

Proposition 1: the svmplectic propertv . 



N 



det L] 
det Li 



deiT{E) = ]J 



(7) 




<N, 
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Proof: it is a consequence of Lq = Ln and of the factorization of T{E) into a 
product of matrices Tk{E), whose inverse is: 

Corollary: if z is an eigenvalue of T{E), then 1/z* is an eigenvalue of T{E*). 
For real E, both z and 1/z* are in the spectrum of T{E). If T{E*) = T{Ey, 
then both z and 1/z arc in the spectrum of T{E). 
These statements summarize in the useful identity: 

det[r(£;) -z]= z^^ detT(O) det[T{E*) - 1/z*]* (9) 



Let u, with blocks Ui, . . . ,Un, be an eigenvector of with eigenvalue E. 
Then, by eq. (3) and after imposing the boundary conditions (4): 




which means that z is an eigenvalue of T{E) with an eigenvector of blocks 
zui and u^. However the converse is true: given an eigenvector of T(E) with 
eigenvalue z one reconstructs via products of matrices Tk{E) the whole eigen- 
vector of TL{z) with eigenvalue E. Therefore det[E — Ti.{z)\ = if and only if 
det[T(_E) — z] = 0; this duality among eigenvalues is made precise as an identity 
among characteristic polynomials [8]: 



Proposition 2: the spectral duality . 

Aet[E - H{z)] = {-^z)'^' Aet[LiL2 ...Ln] dct[r(£;) - z] (11) 

Proof: we must show that det [T(i?) — z] is a polynomial in E of degree NM 
with leading coefBcient (— z)*^det[ii . . .Lm]~^- To this end, we first consider 
the leading terms in E of both sides of cq.(9): the equality implies that the 
leading term of deX\T{E) — z\ is proportional to z^ . Next we derive the following 
leading block stucture of T{E): 



The leading term in E of dot [T(i?) — z/] with the constraint of being proportional 
to z^ , is provided by the diagonal factors det{E^Q) and det(— zl). • 



Corollary: If Imi? 7^ then T{E) has no eigenvalues on the miit circle. 
Proof: for |z| = 1 the matrix H(z) is Hermitian, and for ImE ^ it is always 
det[£; - W(z)] 7^ 0. By duality, this implies det[T{E) - z]^0. • 



Notes: in [9] I provided a representation of T{E) — z in terms of the corner 
blocks of the resolvent of 'H{z). I also state the spectral duality for the matrix 
TtT. 
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In the tridiagonal case (M = 1), blocks are numbers. If A = Li...Ln the 
spectral duality simplifies to the known expression 

det[E-n{z)]= XtiT{E)- Xz- X*- (12) 



3. The spectrum of exponents 

Let us fix E real or complex and denote as Za{E) = 6?"+*'''", o, = 1 . . . 2M, the 
eigenvalues oiT{E). The real numbers E,a(.E) are the exponents of the transfer 
matrix. For real E the symplectic property (8) assures that exponents come in 

pairs ±ia{E). 

The property | deiT{E)\ = 1, see (7), implies: 

2M 

Y.UE)=^ (13) 

o=l 

When considering an ensemble of random Hamiltonian matrices with tridiagonal 
block structure, one is often interested in the corresponding ensemble of transfer 
matrices. Being a product of N random matrices, a transfer matrix develops 
exponents £,a{E) that asymptotically grow linearly in the length N [20], with a 
coefficient known as the Lyapunov exponent: 

7a(i?) = lim hum (14) 

N—^oo iV 

For tridiagonal random matrices there is just one pair of opposite Lyapunov 
exponents, that can be evaluated with the Herbert- Jones-Thouless formula [21], 
with the knowledge of the average density of eigenvalues: 



j{E) = const + j dE'p{E') log \E - E'\ 



The extension to a complex value E is discussed in [22]. 

For the Anderson model [23] or Band Random Matrices [24] there are several 
numerical studies of Lyapunov spectra. In these cases of great physical interest, 
transfer matrices are derived from the Hamiltonians and the analytic approach 
is difficult. 

It is interesting that an analytic formula relating the distribution of the ex- 
ponents to the spectrum of the Hamiltonian is possible. Note that the following 
statements are true for a single and general block-tridiagonal Hamiltonian. 
We shall deduce several results from 



Proposition 3: 

■ log I Aet[E - H(e«+'^)]1 = ^log | detL,| + i |^ - UE)\ (15) 



_ „ , 2 

i=\ a=l 
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Proof: The duality relation (11) gives: 

log|det[£;- W(e«+*'^)]| -^log|detLi| = 

i 

2M 

2 

2M 

2 



= -M^ + i Vlog|e«"+^'^» _e«+»¥'|2 

I lii^a + l0g[2c0sll(^„ - - 2cOs((pa " ^)]} 

2M cx> 2M „/ s 

4ei^--^i-ee ^" ^"r^^ ^-^'^-^' (16) 

a=l ^=1 o=l 

We used eq.(13) and the formula (see eq. 1.448.2 in [25]): 

log[2cosha;-2cosy] = |x| - 2 ^ ^^e"^'^! 



Eq. (16) is the Fourier expansion of log | dct[i? — H(e^+*''')] |, which is a periodic 
function of (p. The constant mode is just the proposition. • 

In the special case ^ = 0, the matrix TL{e'^'^) is Hermitian, and (15) yields a 
formula for the sum of positive exponents {E can be complex): 
Proposition 4: 

^$„(i?) = -^log|detL,|+ /■'"|^log|det(£;-W(e'^))| (17) 

By taking the derivative in ^ of (15) we obtain the spectral counting function, 
which counts the exponents less than ^, for any complex value E\ 

2M 

AA(e.(i?)<0 = E^(^-^'^(^))= (18) 

a=l 

M+^£^glog|det[E-W(e«+^^)]| (19) 

We now write | det[i? — 'H{z}] \ in terms of the eigenvalues En{z) and their com- 
plex conjugate, which equal En{l/z*), and evaluate the derivative: 

^log|det[£;-H(e«+'n]l = 
1^/ 1 dEn{z) , 1 dE„{zy 



El 1 i^^ny^) 

. ^ \E-Er,{z) ^ E*-E„{z)* 

1^/ 1 dEr^jz) 1 dE„{zr \ 

2iZ^\E - En{z) dip E*-En{z)* dip J ^ ' 
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Let us denote with N-i.{E) and N-{E) the numbers of positive and negative 
exponents of T{E). We have: 

Proposition 5: N+{E) = N_{E) 

Proof: for real E we know that the exponents come in pairs ±^a{E), because 
of proposition 1. Let us consider the case ImE ^ 0. 

As a consequence of duahty we derived that no eigenvalue of T{E) is on the 
unit circle, thus all exponents are non-zero: N- + N+ = 2M. Therefore, if 
wc set ^ = in (19), the left term is 7V_. Wo now show that = M or, 
equivalently, that the integral in (19) vanishes for ^ = 0. In the expression (20) 
the eigenvalues £'„(e'*') are real periodic functions of ip in [0, 27r], then 



As a function of (p (^ is fixed) an eigenvalue En{e^^'^'^) makes a loop 7„ in the 
complex E plane. The loop 7* of E* is specular with respect respect of the real 
axis. Integration in (fi of (20) yields Cauchy integrals 



The first integral is the winding number of the (oriented) trajectory 7„ of 
En{e^^^''^) around the value E. The second integral is the winding number 
of 7* around E* and has opposite sign because of opposite orientation. 
We then obtain a nice geometric result: 

Proposition 6: M{£,a{E) < ^) = M + W{E) 

The number of exponents of T{E) less than ^ equals M plus the total winding 
number of loops of eigenvalues £'„(e^+"''), —it < if < n, that encircle E 

As I said, these formulae hold for a single general block tridiagonal Hamilto- 
nian matrix and its transfer matrix. For a statistical ensemble of Hamiltonians 
one performs ensemble averages in place of a phase average, and deals with Lya- 
punov exponents. Souillard [20] obtained the following formula for the positive 
Lyapunov exponents, which is the statistical analogue of (17): 



p{E) is the ensemble averaged spectral density of the Hermitian Hamiltonians. 
I am not aware of any statistical analogue of (19). 

4. Bands, arcs and energy level motion. 






(21) 
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Spectral duality provides information on the positions of the eigenvalues of 7i (z) 
and their motion under variations of the boundary parameter z. 
It is useful to introduce the notion of discriminant. If Za{E) is an eigen- 
value of T{E), the discriminant is the eigenvalue of T{E) + T{E)~^ with same 
eigenvector: 

Aa{E) = Za{E) -\ i— = 2 cosh COS fa + ^i smh simpa (22) 

Za{E) 

Since the symplectic property implies that 1/z* is an eigenvalue of T{E*), it is, 
in general: 

A,{E*) = A^iE)* (23) 



Let us begin with the simpler case where Hi, Li are real matrices. Then 
T{E*) = T{E)* and the eigenvalues of the transfer matrix come in pairs Za and 
a = 1 • • • M . Moreover, if E is real the characteristic polynomial of T(E) 
has real coefficients and roots come also in pairs Za, z*. 
The spectral duality reads: 



1\ f 1 



n 



ZaJ \ Z 



det[Lx...Ln]~^ det[E -n{z)] (24) 



Therefore the M equations 



A,{E)=z + - (25) 
z 



provide the NAd eigenvalues of TC{z), which arc naturally classified in subsets 
with label a. We consider two cases: \z\ = 1 and z real. 

When z = e'*' all eigenvalues of H(e"^) are real periodic functions of — tt < 
(f < TT. Each equation Aa{E) = 2 cos 95 provides a number Ua > of real 
solutions, and ni + . . . + um = NM. This means that y = Aa{E), as a function 
of the real variable E, crosses Ua times the strip — 2 < y < 2 parallel to the E 
axis. No extrema are allowed in the strip, for all branches must cut Ua times the 
lines y = 2cosyj, to ensure that 7i(e*'^) has NM real eigenvalues for all cp. All 
branches of the functions Aa{E) cross their bands in the eigenvalues of H{i). 
See Fig. 1. 

By changing ip the eigenvalues span bands in the real axis 

E = Ea,j{'f) a= 1...M, J = l...n„ 

with extrema given by the eigenvalues of the periodic {^p = 0) and antiperiodic 
(v? = tt) Hamiltonians. The velociy of level motion is 

dE _ 2siny 
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In the turning points the velocity vanishes. The second derivative is known as 
the curvature [1, 26, 27]. This dynamics is invariant under the "time-reversal" 
operation ip —ip, corresponding to the transposition of the Hamiltonian ma- 
trix. 

Bands with same label a may at most share an extremum, while bands related 

to different a may overlap. 

When a branch of and a branch of A^' cross inside the strip in a point 
{E,2cosip), there is a crossing of two eigenvalues of H{e^'^) (and a collision of 
two pairs of eigenvalues of T{E) in the unit circle). This is a highly non-generic 
occurrence for one-parameter Hermitian matrices. 

Energy bands of Hermitian periodic tridiagonal matrices (M = 1) were studied 
in [28] 



Fig.l Energy bands of a pentadiagonal matrix (M=2), solution to Aa{E) = 
2 cosy, a = 1, 2, —IT < ip < n. 

When z = e^^, with g > 0, the matrix T-l{z) is real. The equations Aa(E) = 
2cosh.{Ng) provide all NM eigenvalues in subsets labelled by a. They imply 
the equation ^a{E) = ±Ng. For finite g and large N, given that the exponents 
grow linearly in N with a coefficient that defines the Lyapunov exponent, the 
eigenvalues of the non-Hermitian matrix 'H{z) distribute along M level lines 

la{E) = g (27) 

The distribution along arcs is shown in Fig. 2, for large z. For small z the 
pattern of the eigenvalue distribution is complicate and intertwined. 
For g small enough, in continuity with the description given for the periodic case, 
the eigenvalues still belong to the real axis, outside their periodicity bands. By 
increasing g, eigenvalues with the same label a pairwise approach until a pair 
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condenses. Correspondingly, an extremum is reached for the function Aa{E) 
and the pair of eigenvalues acquires opposite imaginary parts. There is a critical 
value ga of g for this to happen for each label a: 



K{E^) = A,{E^) = 2cosh{Ng,) (28) 
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Fig. 2 The 600 eigenvalues of an eptadiagonal matrix (M=3, N=200), with diag- 
onal uniform disorder \ai\ < 0.5, unit off-diagonal elements, and z = 20 (left), 
z = 390 (right). 

Finally, let us mention the case where Hi = hJ and Li are not real matrices. 
Another spectral identity is needed, which follows from spectral duality and the 
symplectic property [8] . It holds for any value of z or E in the complex plane: 

det[T{E)+T{E)-'^ - (z + -)] = \ det Li . . . L^}'^ det[E - n{z)] det[E - n{-)] (29) 

z z 

For z = e**^ the right side is zero in 2NM real solutions, Ei(ip) and Ei{—ip), 
i = 1 . . . NM, that span the same NM bands as (p varies in —n, tt. They are 
degenerate in if — {y — 2) and ip = n {y = —2). 

The strip —2 < j/ < 2 is thus crossed by 2NM branches y = Aa{E) which join 
pairwise at the boundaries y = ±2 of the strip. Each pair, when intersected 
with the line y = 2 cos ip, determines the same band which the eigenvalue E{(fi) 
covers with different speeds in the two directions. 

For z = ±e^^ we have again condensation of pairs of eigenvalues, which no 
longer become complex conjugate pairs. For large N the eigenvalues move into 
M level curves 7o(£^) = g (the case g < leads to a different set of curves since 
in this case Lyapunov exponents need not be opposite pairs). 

Conclusion 

The spectral duality is a simple identity that links the eigenvalues of a matrix 
with block-tridiagonal structure to those of the corresponding transfer matrix. 
I have deduced some analytic properties for the exponents (sum of exponents, 
= N^, counting function). They involve a phase average on eigenvalues of 
the block-matrix. 
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The large N stability of exponents allows to describe the distribution of complex 
eigenvalues of the block matrix along arcs. Discriminants classify real eigenval- 
ues (periodic case) in bands, and govern their motion and collisions. 
These exact properties are hoped to allow a more analytic approach to the 
difficult study of Lyapunov spectra of transfer matrices, that are derived from 
an ensemble of random Hamiltonians. They also extend to block matrices some 
results which were known for purely tridiagonal matrices. 
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